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ABSTRACT 

Context. The slope of the locally measured spectrum of cosmic rays varies from 2.8 for protons with energies below 
200 GeV down to 2.5 for heavy nuclei with energies in the TeV-PeV range. It is not clear if the locally measured slope 
values are representative for those of the overall population of Galactic cosmic rays and if the slope of the cosmic ray 
spectrum varies across the Galaxy, e.g. in response to the variations of the star formation rate. 

Aims. We use the data of Fermi Space 7 -ray Telescope to derive a measurement of the slope of the cosmic ray spectrum 
across the Galactic Disk and to compare it with that of the Large Magellanic Gloud cosmic rays. 

Methods. A special choice of the background estimation regions allows us to single out the neutral pion decay component 
of the 7 -ray flux in the energy range above 10 GeV and to separate it from (a) emission from the local interstellar 
medium around the Solar system and (b) from the inverse Gompton emission produced by cosmic ray electrons. 
Results. The spectrum of the pion decay 7 -ray emission from the Galactic disk in the energy band 10 GeV - 1 TeV has 
the slope ~ 2.4. There is no evidence for the variation of the slope with Galactic longitude / distance from the Galactic 
Gentre. The slope of the spectrum of cosmic rays derived from the 7 -ray data, ~ 2.45, is harder than the slope of the 
locally observed cosmic ray proton spectrum. Pion decay emission from a powerlaw distribution of cosmic rays with 
the same hard slope also provides a fit to the 7 -ray spectrum of the Large Magellanic Gloud. 

Conclusions. Identical and hard slopes of the spectra of cosmic rays in the Milky Way and in the Large Magellanic Gloud 
are consistent with a straightforward theoretical model in which cosmic rays are injected by shock acceleration with the 
spectrum with the slope P ~ 2...2.1 which is subsequently modified by AT = 1/3...1/2 by the energy-dependent escape 
of cosmic rays through the turbulent Galactic magnetic field. Deviations of the locally measured cosmic ray spectrum 
from the average Galactic spectrum are explained by the discrete distribution of the cosmic ray sources in space and 
time. 
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1. Introduction 

Our knowledge of the properties of Galactic cosmic rays 
(CRs) is based mostly on the measurements in a single point 
in the Milky Way galaxy: the position of the Solar system. 
From these local measurements we know that the spectrum 
of CRs is an approximate powerlaw with the slope ranging 
from 2.8 below 200 GeV to 2.6 ... 2.7 in the TeV-PeV energy 
range, ending with the ” knee” feature above the PeV energy 
(|01ive et al.ll20l4l . 

Numerous theoretical models were developed over time 
to explain the spectral properties of the cosmic ray flux 
(jBerezinskii et all 11984 iBlasil l2013ll . An implicit assump¬ 
tion of these models is that the approximate 2.7 slope pow¬ 
erlaw spectrum is ” universal”, in the sense that it is valid all 
over the Milky Way galaxy. The same assumption is com¬ 
monly adopted in the interpretation of the data on 7-ray 
emissi on from Galactic CR population (|Ackermann et aO 
l2012bll and in the modelling of the extragalactic 7-ray back¬ 
ground produce d by a population of star -forming galaxies 
in the Universe (jAckermann et ani2012ari . 

An assumption that the Galactic CR spectrum has the 
slope 2.7 everywhere across the Galaxy is, however, far from 
being obvious. None of the theoretical models which ’’ex¬ 


plains” the spectral slope is developed from the first prin¬ 
ciples. The match between the observed spectral slope and 
the model calculations is achieved via adjustment of phe¬ 
nomenological parameters, e.g. of the slope of the injection 
spectrum of CRs in the interstellar medium, Tinj and the 
modification of the slope by the energy dependent escape 
of the CRs, charac terised by the energy sl ope of the diffu¬ 
sion coefficient, S (iBerezinskii et allll984ll . It is well pos¬ 
sible that the slope of the locally observed CR spectrum 
is not identical to that of the average Galactic CR spec¬ 
trum. Variations of the star formation rate over the last 
tens of millions of years and discreteness of the cosmic ray 
source distribution could lead to fluctuat ions of the shape 
of the CR spectrum in space and time ( lAharon iaiil |2004 


iNeronov fc Semik^l2012l: iKachelriess et al.ll2015ll. 


The only possibility for the measurement of the spec¬ 
trum of CRs at differe nt locations in the Galaxy is via the 
use of the 7-raY data (fAharonianI 12004 lAckermann et all 
l2012bl : INeronov et al.l[^12ll . CRs interacting in the inter¬ 
stellar medium (ISM) produce neutral pions which decay 
into photons. The spectrum of the pion decay 7-rays has 
the slope which is nearly identical to that of the parent CR 
spectrum. 
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An immediate problem for the measurement of the 
CR spectrum from the 7 -ray data is that the 7 -ray 
flux from the Galaxy contains, apart from the pion de¬ 
cay flux, also the flux produced by the CR electrons 
via the Bremsstrahlung a nd inverse Compton mechanisms 
(|Ackermann et ani2012bll . As a consequence, derivation of 
the CR spectrum from the 7 -ray data should include a de¬ 
tailed modelling of the three contributions, or a method of 
separation of the pion decay contribution to the flux from 
the Bremsstrahlung and inverse Compton contributions. 

In what follows we develop a method which allows to 
separate the pion decay emission from the Bremsstrahlung 
and inverse Compton emission in the flux coming from the 
distant parts of the Disk of the Milky Way galaxy. Using 
this method, we are able to measure the spectrum of cosmic 
rays at different locations in the Disk, at different distances 
from the Centre of the Galaxy. We show that the average 
Galactic CR spectrum in the TeV energy band is harder 
than the locally observed spectrum of the proton compo¬ 
nent of the CR flux. We also demonstrate that the aver¬ 
age Galactic CR spectrum slope does not vary significantly 
across the Galaxy and it is also identical to the slope of 
CRs filling the disk of the Large Magellanic Cloud (LMC) 
galaxy. This suggests that the measured average slope of 
the CR spectrum is ’’universal”, in the sense that it is de¬ 
termined by the basic physical laws. It is this slope which 
could be explained by the straightforward theoretical mod¬ 
els of of CR injection and propagation. 

2. Data Analysis 

In our work we use 6.7 years of Fermi/LAT data (from 
August, 4th, 2008 to December 2014). We process the data 
using the most recent version of the Fermi Science Tools 
Software v9r33p0 and the P7REP response function^ The 
spectra of the sources are calculated using two comple¬ 
mentary techniques: the binned likelihood analysis and the 
aperture photometry. 

The binned likelihood analysis is based on the fitting 
of a model of diffuse and point source emission within a 
’’region of interest” to the data. The spatial extent of the 
region is 13° radius at energies > 300 MeV and 20° ra¬ 
dius in the 60-300 MeV energy range. The spatial model 
includes diffuse Galactic and extragalactic backgrounds 
and the sources from the 4 year (3 FGL) Fermi catalog 
(|The Fermi-LAT Collaborationll^OlSD . In each energy bin 
we fix the spectral shape of each source to be a power law 
with index —2. The spectral shapes of diffuse backgrounds 
are given by the corresponding templates. The normalisa¬ 
tions of fluxes of all sources and diffuse backgrounds are 
treated as free parameters during the fitting. 

The aperture photometry method of spectral analysis 
uses a selection of 7 -ray events from pre-defined signal and 
background regions with the gtbin / gtmktime tools and 
the calculation of the exposure for the selected regions with 
the gtexposure tool. All the sources have the spatial extent 
much larger than the size of the Point Spread Function 
(PSF) of LAT in the energy range above 1 GeV, where 
the aperture photometry analysis is performed. Taking this 
into account, we use the gtexposure tool with the switch 
appcor=N0, which precludes the correction of exposure for 
the PSF effects. We remove the contribution of bright 

^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 


point sources to both the signal and background regions 
by excluding circles of the radius 0.5° around identified 
puls ars and active galactic nuclei from the Fermi cata¬ 
log (|The Fermi-LAT Collaboration! [^OlhTl and around the 
Galactic Centre source. In the analysis of the signal from 
the Galactic Plane we explicitly do not remove the uniden¬ 
tified sources from the catalog, because a large fraction of 
these sources could be just the inhomogeneities of the the 
overall diffuse emission from the ISM. We also do not re¬ 
move the supernova remnants because they could be consid¬ 
ered as a part of the overall 7 -ray signal from the interacting 
Galactic CRs. 

3. The Galactic Disk 

The Galactic diffuse emission is a sum of the pion de¬ 
cay, Bremsstrahlung and inverse Compton contributions. 
Relative importance of different contributions to the 7 -ray 
flux depends on the photon energy and on the location of 
the emission region in the Galaxy. 

Our goal is to derive a measurement of the spectrum of 
proton / nuclei CRs residing at different locations in the 
Galactic Disk from the spectrum of the pion decay emis¬ 
sion from the Disk. There are two obvious obstacles for 
such a measurement. First, the pion decay flux has to be 
separated from the Bremsstrahlung and inverse Compton 
fluxes. Next, the flux of the pion decay emission produced 
at large, several kiloparsecs distances from the Sun has 
to be separated from the flux of pion decay and inverse 
Compton / Bremsstrahlung emission form the local inter¬ 
stellar medium. 

The spectrum of the Bremsstrahlung is soft and its con¬ 
tribution to the 7 -ray flux is negligible in the energy band 
above 10 GeV (jAckermann et al.l l2ni2bll . Thus, reducing 
the energy range of the analysis to > 10 GeV should pro¬ 
vide a signal which is largely free from the Bremsstrahlung 
contribution. 

The spectrum of inverse Compton emission is harder 
and the inverse Compton flux could be significant up to 
the highest energies accessible for the Fermi/LAT. The in¬ 
verse Compton flux is strongest close to the direction of the 
Galactic Plane and toward the inner Galaxy, because of the 
higher density of the interstellar radiation field in the inner 
part of the Galactic Disk. 

Fortunately, the pion decay 7 -ray emission is still more 
concentrated toward the Galactic Plane because its volume 
luminosity scales with the density of the interstellar gas. 
The luminosity is falling exponentially with the increase of 
the altitude above / below the Galactic Plane because the 
density of the Galactic Disk is falling exponentially with 
the scale height about 100 pc. 

This difference leads to a significant difference of 
the Galactic latitude pr ofiles of the pioii decay and in¬ 
verse Compton emission (lAckermann et aDl 2012 bll . see Fig. 
dl right panel. To separate the pion decay and inverse 
Compton signals, we take the source signal from a nar¬ 
row strip | 6 | < 1.5° around the Galactic Plane and esti¬ 
mate the background from two narrow strips 2° < | 6 | < 3° 
as it is shown in the left panel of Fig. [T] The pion de¬ 
cay emission profile is sharply peaked toward 6 = 0 °, com¬ 
pared to the inverse Compton profile. Integrating the pion 
decay and inverse Compton fluxes in the | 6 | < 1.5° and 
2° < | 6 | < 3° regions one could find that the choice of the 
source and background regions shown in the left panel of 
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Fig. 2. Spectrum of emission from the |Z| < 90°, |&| < 1.5° 
part of the Galactic Plane. Thick / thin errorbars show 
the statistical / systematic error. Grey thick line shows the 
best-fit powerlaw with the slope r.y = 2.42 in the energy 
range above 10 GeV. Dashed line shows the spectrum of the 
neutral pion decay emission produced a powerlaw proton 
spectrum with the slope Pp = 2.45. 

Fig. [T] rejects cs 80% of the inverse Gompton signal within 
|6| < 1.5° strip, while leaving > 50% of the pion decay flux. 
In this way the pi o n deca y flux, which (within the model of 
lAckermann et all (|2012blD was already an order of magni¬ 
tude higher than the inverse Compton flux within the signal 
extraction strip, gets boosted by another half-an-order. 

The same choice of the signal and background regions 
allows also the suppression of the flux from the local ISM. 
The angular extent of the source occupying a 100 pc wide 
layer around the Galactic Plane is about 10° at the distance 
1 kpc from the Sun and is just about 3° for the flux coming 
from the distances beyond 3 kpc from the Sun position. 
The highest star formation rate in the Milky Way is in 
the innermost part of the Galactic Disk at the distances 
> 4 kpc from the Sun. This innermost part of the Galactic 
Disk spans a strip |6| < 3° on the sky. Choosing both the 
signal and background regions within a region where the 
flux from the local ISM does not vary removes the local 
ISM contribution. 

The spectrum of diffuse emission from |6| < 1.5°, |/| < 
90° strip extracted using the aperture photometry method 
is shown in Fig.j^J The spectrum above 10 GeV could be fit 
by a powerlaw with the slope P.^ = 2.42±0.03stat±0.12sys(. 
In this energy band the flux is dominated by the pion decay 
emission. The spectrum could be fitted with a model pion 
decay spectrum calculated using the code of iKamae et al.l 
(|2006r) for the powerlaw proton spectrum with the slope 
Pp = 2.45. 

Below 10 GeV the y-ray spectrum gets softer. This 
softening is expected at least due to a non-negligible con¬ 
tribution of the electron Bremsstrahlung to the photon 
flux. Alternatively, the pion decay spectrum itself could 
soften below 10 GeV is the CR spectrum gets softer below 

150 GeV, as does the local CR spectrum. Our analysis 
method does not allow to distinguish between these two 
possibilities. 

It is not clear a-priori if the average slope of the CR 
spectrum depends on th e distance from the c entre of the 
Galaxy as suggested by iGaggero et al.l (l2015ll or it has a 
’’universal” slope. It is also not clear if the slope depends 
on the relevant physical parameters like e.g. the star forma¬ 
tion rate. Fig. [3] shows the slope of the y-ray spectrum of 
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Fig. 3. Galactic longitude profile of the flux (bottom panel) 
and slope (top panel) of y-ray emission from the Galactic 
Disk in the energy band above 10 GeV. Horizontal black 
line shows the best fit spectral slope. Orange Ithin line and 
the hatched region show the slope of the y-ray spectrum of 
the LMC. Green and olive data points at I ~ 80° show the 
slope of the spectrum of Cygnus region extracted using the 
likelihood and aperture photometry methods. Light blue 
data points show the slopes of the spectra of unidentified 
sources from the HESS Galactic Plane survey. 

the |6| < 1.5° strip as a function of Galactic longitude. One 
could see that there is no evidence for the dependence of 
the slope on the Galactic longitude or, equivalently, on the 
distance from the Galactic Centre. There is also no correla¬ 
tion between the slope and the overall flux above 10 GeV, 
as one could see from comparison of the upper and lower 
panels of Fig. |31 

3.1. Cygnus region 

The method of calculation of the pion decay contribution 
to the diffuse flux from the Galactic Disk developed in the 
previous section does not work in the outer Galaxy, in the 
Galactic longitude range |Z| > 90°. This is because the con¬ 
tribution of the distant parts of the Galactic Disk at several 
kiloparsecs distances becomes much smaller than the flux 
from the local ISM. In this case, estimation of the flux from 
a narrow strip around 6 = 0° plane catches the degree-scale 
fluctuations of the flux from the local Galaxy, rather than 
the flux from the distant parts of the Galactic Disk. 

A measurement of the spectrum of CRs in a distant part 
of the Galactic Disk is still possible in one particular case, 
namely in the direction of Cygnus region. This direction is 
tangent to the Galactic arm passing close to the Sun. An 
enhancement of the diffuse emission from the direction of 
Cygnus region is commonly attributed to a superposition 
of the projected y-ray emission from the local Galactic arm 
and the flux produced by the nearby active star formation 
region at ~ 1.5 kpc distance in this direction. Emission 
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b, degrees 


Fig. 1. Left: Fermi/LAT count map of the inner Galaxy region above 1 GeV. White solid boxes show the 3° x 20° regions 
of signal extraction, dashed green boxes show the regions of background estimate. Right: Galactic lat i tude pr ofiles of 
the 7 -ray emission from pp interactions and of the inverse Gompton emission, from lAckermann et ;iI1 (l2012bfl . Dense 
and light shading show the ranges of Galactic latitudes from which the signal and background are extracted for the 
calculation of the spectrum. 



Fig. 4. Fermi/LAT count map of Gygnus region in the en¬ 
ergy band above 1 GeV. White circle shows the extent of 
the Gaussian model template for the spectral extraction. 
White solid and green dashed boxes show the regions of 
signal and background extraction for the calculation of the 
spectrum using the aperture photometry. 



E, GeV 

Fig. 5. Spectrum of the Gygnus region extracted from a 2° 
Gaussian spatial template (thick data points) and from an 
5° X 5° box. Grey dashed line shows a fit with a model of 
pion decay emission from the powerlaw CR spectrum with 
the slope Tcr = 2.45. Solid grey line shows a powerlaw fit 
to the box spectrum. 


from this nearby star forming region is known to have a 
hard spectrum with the slope close to L.^ ~ 2.2, suppos- 
edly due to a recent in jection of freshly accelerated CRs 
(|Ackermann et ^l201lH . 

The extended emission source in the Gygnus region, re¬ 
ported in the Fermi catalog, is described by a Gaussian spa¬ 
tial template with the width 2°, shown in Fig. |4l Spectrum 
of the source with such spatial morphology is shown in Fig. 
O The spectrum above 3 GeV is well described by a power- 
law with the slope T.,, = 2.47T0.09. This slope is softer than 
the slope reported bv lAckermann et al. (120111). because we 
have used, contrary to lAckermann et al.l ~( 201l!l . the stan¬ 
dard Galactic diffuse emission background for the spectral 
modelling. In such standard model the extended source in 
the direction of Gygnus includes the entire emission from 
the Galactic arm behind the nearby star formation region, 
rather than only a contribution from the nearby star for¬ 
mation region at 1.5 kpc distance. 

As a cross-check, we extract the spectrum of the Gygnus 
region using the aperture photometry method and taking 
the source region to be a square box of the size 5° x 5° 
centered at I = 80°, & = 1°, see Fig. 01 Similarly to the 
Galactic Plane spectrum case, we suppress the contribu¬ 
tion to the flux from the local interstellar medium and from 
the inverse Compton emission by an appropriate choice of 


the background estimate regions. Both the local interstel¬ 
lar medium and inverse Compton contribution do not vary 
significantly along the Galactic Plane in the Galactic longi¬ 
tude range 70° < I < 90°. Choosing the boxes of the same 
size as the signal box, but displaced toward smaller and 
larger 1, as it is shown in Fig. 0] suppresses the two com¬ 
ponents of the diffuse emission, leaving the pion decay and 
Bremsstrahlung flux from the Gygnus region unaffected. 

The spectrum of emission from the 5° x 5° box is shown 
in Fig. by the thin olive colour data points. The slope 
of the spectrum is r.y = 2.42 ± 0.06, which is consistent 
with the slope of the spectrum extracted using the Gaussian 
template. 

The emission from the Gygnus region is dominated by 
the pion decay flux, as it is clear from the presence of the 
low-energy cut-off in the sp ectrum at the pion prod uction 
threshold ~ 100 MeV (lAckermann et all l2013ll . The 
slope of the pion decay spectrum in the Gygnus region 
is identical to the slope of the spectrum elsewhere in the 
Galactic Plane. This confirms that the slope P.^ ~ 2.42 is 
characteristic for the slope of the CR spectrum produced 
by star formation activity and that it is largely indepen¬ 
dent of the parameters of the interstellar medium. Indeed, 
the properties of the interstellar medium in Gygnus region 
and in the local Galactic arm behind the Gygnus region 
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Fig. 6. Fermi/LAT count map of LMC field in the energy 
band above 1 GeV (left) compared to the Ha map of LMC 
(right). Smaller white circle in the left panel shows the ex¬ 
tent of the Fermi/LAT Gaussian template for the spatial 
model of LMC. Larger circle shows the extent of the source 
signal extraction region for the aperture photometry calcu¬ 
lation of the spectrum. 



Fig. 7. Spectra of the LMC extracted using the likelihood 
analysis (thick data points) and the aperture photometry 
method (grey thick / thin data points showing the statis¬ 
tical / systematic errors). Dashed line shows a model fit 
with the pion decay spectrum produced by the parent pro¬ 
ton population with the powerlaw spectrum with the slope 
Tcfl = 2.45. Solid line shows a powerlaw fit to the aperture 
photometry spectrum. 


are quite different from those in the innermost arms and in 
the Galactic Bar. However, the slopes of the CR spectra in 
these regions are the same. 


4. Large Magellanic Cloud 


The conjecture on the universality of the spectrum of CRs 
produced by the star formation could be tested via obser¬ 
vations of other galaxies. The closest to us disk galaxy is 
the Large Magellanic C loud (LMC) at the distance 50 kpc 
(|Pietrzviiski et ani2013ll . 

The Fermi/LAT count map of the LMC field is shown in 
the left panel of Fig. [6] The brightest emission is detected 
from the direction of 30 Doradus star forming region, which 
is also visible as a bright excess on th e Ha map of the LMC 
shown in the right panel of Fig. [6] (ICaustad et al.l[l999h . 
The overall extent of the LMC disk is at least 3°, as it is 
clear from the Ha map. However, the y-ray brightest part of 
the LMC is more compact, fitting within a circle of ~ 1.9° 
in radius. Apart form the 30 Doradus region, significant 
flux comes from the Constellation HI star forming region 
North of the 30 Doradus. 

The overall spati al model of the ^ - ray e mission from 
LMC, considered by lAbdo. A.A. et al. (|201flll is a sum of 


two gaussians of 1.9° and ^ 0.2° widths. The smaller 
size Gaussian is centred at the position of the 30 Doradus 
region. The spectrum of LMC extracted using this template 
is shown in Fig. [71 Above 3 GeV, the LMC spectrum is 
well described with a power law model with the index T = 
-2.41 ±0.15. 

Similarly to the case of the Cygnus region, we cross¬ 
check the result using the aperture photometry spectral ex¬ 
traction method. We take the signal region to be a circle of 
the radius 3°, centered at RA=79.61°, DEC=-68.53°, which 
corresponds to the position and extent of the Ha disk of the 
LMC, see Fig.|6l We estimate the background from a circle 
of equal radius centered RA=70.82°, DEC=-76.75°, situ¬ 
ated at the same Galactic latitude as the source circle. The 
spectrum extracted using the aperture photometry method 
is shown by thin data points in Fig. [T] It could be fit by 
a powerlaw with the slope T.,, = 2.39 ± O.OQstat ± O.OTgyst, 
which is consistent with the slope of the spectrum extracted 
using the likelihood analysis. 

The overall spectrum in the broad energy range from 
60 MeV up to 100 GeV is well fit with a single pion de¬ 
cay component spectrum with the slope of the parent CR 
spectrum Tc/j = 2.45, shown by the dashed line in Fig. 
[T] Contribution of the electron Bremsstrahlung component 
to the spectrum is constrained by the upper limit in the 
60-100 MeV band, where the Bremsstrahlung flux should 
dominate over the pion decay flux. Similarly to the Galactic 
Plane and Cygnus region spectra, the Bremsstrahlung flux 
in the energy range above several GeV is negligibly small. 

The slope of the y-ray spectrum of the LMC above 
3 GeV is consistent with the slope of the spectrum of y-ray 
emission from the Milky Way disk, see Fig. [31 This sug¬ 
gests that the spectrum of emission from the LMC disk is 
dominated by the pion decay flux from the proton / nuclei 
CR interactions in the LMC disk, similarly to the Milky 
Way emission in the same en ergy band. This concl usion 
is supported by the analysis of [Foreman et al.l (l2015ll who 
find that the pion decay flux dominates the y-ray emission 
from the LMC except possibly for the 30 Doradus region 
where electron Bremsstrahlung and inverse Compton emis¬ 
sion could compete with the pion decay emission. 


5. Other star forming galaxies 

Apart from the Milky Way and LMC, several other nor¬ 
mal star forming Galaxies are detected by Fermi/LAT. The 
conjecture of the universal hard slope of the CR spectrum 
produced by the star formation should, in principle, hold 
also for these galaxies. 

We have extracted the spectra of the Small Magellanic 
Cloud and M31 galaxies and checked that the quality of 
the spectra is not yet sufficient for the measurement of the 
details of the spectral properties of the CR populations in 
these galaxies. 


6. Discussion 

Our analysis of the spectral properties of y-ray signal from 
a 3° wide strip along the Galactic Plane has enabled a 
separation of the pion decay emission component of the 
Galactic diffuse emission from the Bremsstrahlung and 
inverse Compton components in the energy band above 
10 GeV. We have used the measurement of the slope of the 
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pion decay spectrum to measure the slope of the spectrum 
of CRs with energies above 100 GeV across the Galactic 
Disk. 

The slope of the Galactic cosmic ray spectrum turns out 
to be harder than that of the locally measured CR spectrum 
in the same energy band. The slope of the Galactic CR 
spectrum does not exhibit significant variations from the 
inner to outer Galaxy. This points to a ’’universal” nature 
of the hard slope, which is largely independent of the local 
physical conditions in different parts of the Galaxy. 

The same hard slope of the CR spectrum is found also 
in the LMC galaxy. This further supports the conjecture 
about the ’’universality” of the hard slope of the spectra of 
CRs resulting from star formation in galaxies. 

The slope of the average CR spectrum derived from 
our analysis of the Galactic Plane and of the LMC coin¬ 
cides with the slope of the neutrino spectrum in the energy 
band above 10 TeV, reported by the IceCube collaboration 
(jAartsen et al.l l2015ll . This suggests that the astrophysi- 
cal neutrino flux has a contribution from the Milky Way 
(|Neronov et al.l[20l4lNeronov fc Semikod[2014ll . A smooth 
match between the diffuse Galactic 7 -ray emission spec¬ 
trum in the sub-TeV energy band with the IceCube neu¬ 
trino spectrum extending up to the PeV energy range indi¬ 
cate that the hard powerlaw CR spectrum with the slope 
PcR — 2.45 extends over several decades in energy, from 
< 1 TeV up to ~ 10 PeV. 

In the 1-10 TeV energy range, not covered by the 
Fermi/LAT and IceCube data, the data of the ground- 
based 7 -ray telescopes could provide a measurement of the 
characteristic slope of the Galactic CR spectrum. Such a 
measurement is more challenging because of the lack of 
the wide-held telescopes in this energy band. However, an 
assessment of the average CR spectrum in the Galactic 
Disk is still possible using the narrow held-of-view tele¬ 
scopes like e.g. HESS. This is illustrated in Fig. |3l The 
thin light blue data points in the upper panel of this hg- 
ure show the measurements of the slopes of the spectra of 
unidentihed extended sources in t he Galactic Plane , from 
the HESS Galactic Plane survey (|Aharonian et ahl 120051 
120001 I 2 OO 8 II . Fairly tight correlation of the positions of 
the extended HESS survey sources with the still more ex- 
tende d sources detected above 1 00 GeV in the Fermi/LAT 
data (iNeronov fc Semi]^l2012ll indicates that the HESS 
unidentihed sources are parts of larger excesses of the dif¬ 
fuse emission occurring at the positions of over-densities 
in the interstellar medium and/or close to the most recent 
points of injection of CRs in the ISM. If this is so, a spread 
in the slopes of the 7 -ray spectra of these sources is ex¬ 
pected because the dense clumps of the ISM could be situ¬ 
ated closer or further from the recent points of injection of 
CRs (|Aharonianll2004lKachelriess et al.ll2015fl . At the same 
time, the average-over-the-sources slope of the CR / 7 -ray 
spectra should be close to the average slope of the spectrum 
of Galactic CRs. This is indeed the case. The average slope 
of the spectra of sources from the HESS survey is remark¬ 
ably close to the slope of the spectrum of the |&| < 1.5° 
strip, Thess — 2.43. 

The value Fc/j ~ 2.45 of the universal slope of the 
spectrum of CRs produced by the star formation process 
could be understood within a straightforward theoretical 
model, without adjustment of phenomenological parame¬ 
ters. The most commonly considered process of the CR 
production is the shock acceleration process. This pro¬ 


cess results in the slope of the injection spectrum of CRs 
c:: 2...2.1 in the most simple se ttings (iBerezinsk i et all 


1984 : Blandford fc Ostrik^ 1978t Berezhko et al.l 1994 : 


Maikov fc Drury 20011 : iBell fc LucekI 20011 :" Blasi l2013ll . 

Propagation of CRs through the turbulent Galactic mag¬ 
netic field with the subsequent escape of CRs from the 
Galaxy leads to the softening of the CR spectrum from 
Tinj to Fcjj = Tinj + S where 6 is related to the slope 
of the turbulence power spectrum. It is 5 = 1/3 for the 
Kolmogorov turbulence and 6 = 1/2 for the Iroshnikov- 
Kraichnan turbulence. Measurement Fca — 2.45 suggests 
that the injection spectrum of CRs in the ISM has the slope 
Tinj = 2.1, assuming the Kolmogor ov spectrum of the ISM 
turbulence ([Armstrong et al.lll995ll . 

Also straightforward is the understanding of the differ¬ 
ence of the slopes of the average Galactic CR spectrum 
and the CR spectra measurable at fixed locations in the 
Galactic Disk, like e.g. the measurements at the position 
of the Solar system. The spectra of CRs measured in any 
single point typically deviate from the average CR spec¬ 
trum because they are influenced by the discreteness of the 
CR source distribution in space and time. A single source 
which injects about 10^° erg in CRs could significantly alter 
the overall energy density of CRs in the local ISM within 
the region of the volume more than ( 0.1 kpc)^ on the time 
scales of millions of yea rs (|Aharonianll2004lAhn et ^120101 : 
iKachelriess et al]l2015ll . In a similar way, a star formation 
activity episode resulting in formation of a superbubble 
with an energy output 10^^ erg would change the CR spec¬ 
trum in a kpc-scale region a round it for tens of millions 
of years (INeronov et al.ll20ia . In general, the CR spectra 
measurable at locations close to recently active sources are 
harder than the CR spectra in a places where no recent in¬ 
jections of CRs happened on the time span co mparable to 
the e s cape time of CRs from the Galactic Disk ( Aharoni^ 
|2004 INeronov fc Semik^ l 2012 t IKachelriess et al.1 201511 . 
The slopes of the CR spectra shown in Figs.[2j[3j[5l[2]are the 
same only because they are derived from the 7 -ray signal 
collected from sufficiently large regions of the ISM, where 
the assumptions on the continuity and time independence 
of the source distribution are valid. 


7. Summary 

We have used the 7 -ray data to derive a measurement of 
the characteristic spectrum of CRs produced by the star 
formation activity in the Milky Way (Fig. [3]) and in the 
LMC (Fig. ini) galaxies. The slope of the average spectra of 
CRs in these two galaxies, Fca — 2.45, is harder than the 
locally measured spectrum of CRs penetrating in the Solar 
system. Consistency of the measurement of the slopes of the 
CR spectra at different locations in the Milky Way Galactic 
Disk and in the LMC suggests that the spectrum of CRs 
produced by the star formation activity is universal, i.e. it 
does not depend on the details of the physical parameters 
of the star formation process and of the ISM. The value of 
the universal slope of the CR spectrum is consistent with 
the straightforward theoretical model of CR injection by 
the shock acceleration process, which produces the spec¬ 
trum with the slope Tinj — 2 ... 2 . 1 , followed by the energy- 
dependent escape of CRs through the turbulent Galactic 
magnetic field with the ISM turbulence power spectrum. 
Fluctuations of the slopes of the CR spectra measurable at 
fixed points in the ISM are attributed to the discreteness 
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of the CR source distributions and variations of the star 
formation rate in time. 
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